Crossing the Phantom Divide with Parameterized Post-Priedmann Dark Energy 
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Dark energy models with a single scalar field cannot cross the equation of state divide set by 
a cosmological constant. More general models that allow crossing require additional degrees of 
freedom to ensure gravitational stability. We show that a parameterized post-Friedmann description 
of cosmic acceleration provides a simple but accurate description of multiple scalar field crossing 
models. Moreover the prescription provides a well controlled approximation for a wide range of 
"smooth" dark energy models. It conserves energy and momentum and is exact in the metric 
evolution on scales well above and below the transition scale to relative smoothness. Standard 
linear perturbation tools have been altered to include this description and made publicly available 
for studies of the dark energy involving cosmological structure out to the horizon scale. 



I. INTRODUCTION 

Observational constraints on the acceleration of the 
expansion have continued to close in on a dark energy 
equation of state of a cosmological constant w e = — 1 
that delineates the phantom divide. Testing the small 
deviations from that value in the future requires a de- 
scription of the dark energy that allows the equation of 
state to evolve across the phantom divide possibly mul- 
tiple times. 

It is well known that single scalar fields are gravita- 
tionally unstable to such a crossing of the phantom di- 
vide 0, 0, S]- Dark energy that is minimally coupled 
to the matter requires additional degrees of freedom to 
cross the divide stably. While specific models with multi- 
ple fields can be constructed @, B| they are cumbersome 
or impossible to implement in a general analysis of the 
dark energy. 

The usual approach in the literature for finessing such 
cases is to artificially turn off the dark energy pertur- 
bations explicitly or implicitly by limiting the range of 
observables. Doing so violates energy-momentum con- 
servation whenever w e ^= — 1 and leads to inconsisten- 
cies between the Einstein equations for the evolution of 
the metric due to the Bianchi identities which can persist 
even on small scales. Though the impact of perturbations 
tend to be small near w e = — 1, cosmological constraints 
often require the exploration of a large swath of param- 
eter space around the maximum likelihood. Excising the 
instability around the transition provides another, albeit 
rather ad hoc approach [f|. 

In this Brief Report, we show that the so-called pa- 
rameterized post-Friedmann (PPF) approach to describ- 
ing linear metric evolution in a Friedmann-Robertson- 
Walker (FRW) universe provides a simple solution to this 
dilemma. The PPF framework was ostensibly introduced 
for describing modified gravity theories under a metric 
framework with strict local conservation of energy and 
momentum Q. As such it also applies to dark energy 
models Q and in particular the class of models which 



have a well-defined Jeans scale under which the dark en- 
ergy is smooth compared to the dark matter. This frame- 
work also has the benefit of being an exact description for 
the metric evolution well above and well below this scale 
and hence provides a very well-controlled approximation 
that is simple to implement in an Einstein-Boltzmann 
linear perturbation code. 



II. PHANTOM DIVIDE AND SCALAR FIELDS 

Minimally coupled scalar field dark energy models that 
evolve across the phantom divide require new internal 
degrees of freedom to maintain gravitational stability. To 
see this fact, consider the conservation equation for the 
momentum density (p e u e )i = T° i (sec, e.g. @ for an 
explicit derivation) 

^Pe 

u' = (3w e - l)u e + kn — - + (1 + w e )k H A, (1) 

Pe 

where ' = d/dhia, w e = p e /p e , ku = k/aH, and A is 
the gravitational potential in an arbitrary gauge. 

The relationship between the pressure and density fluc- 
tuation defines a sound speed. For a single scalar field 
with kinetic and potential degrees of freedom, this re- 
lationship is most simply described in a coordinate sys- 
tem that comoves with the dark energy such that the 
momentum density and transverse spatial metric fluctu- 
ations vanish]^, 0]. From an arbitrary gauge, this quan- 
tity is obtained by a gauge transformation that changes 
the time slicing 



<5p (rcst) = 6p t 
which defines a sound speed 
2 _ S P 



Sp e + 3p e —^. 

k H 



oPe U e 
Pe k H 



(rest) 



5p> 



(rest) 



(2) 



(3) 



2 



bringing the momentum conservation equation to 



W e + C 2 S - ^ - i j U e + kHcj.Se + (1 + W e )kH-A , 



where <5 e = 5p e /p e - 

For a single scalar field, the rest or zero momentum 
gauge corresponds to time slicing where the field, and 
hence the potential energy, is constant leaving the en- 
ergy density and pressure to be defined by fluctuations 
in the kinetic energy. For a canonical kinetic term c 2 = 1 
representing the familiar kinetic energy dominated equa- 
tion of state of such scalars. 

The dark energy system is completed by the continuity 
equation 



S' e + 3{C 2 S - We)5 e + 9 



Pe 



k'H 



-k H u e - (1 + w e ){k H B + 3H' L 



(4) 



where B is the space-time piece and Hl the space-space 
curvature piece of the metric fluctuations in an arbitrary 
gauge [13 • 

Taking c 2 s > makes dark energy perturbations Jeans 
stable in the regime knc s 3> 1 so long as p' e / 1 p' e remains 
finite. In the matter dominated epoch, matter density 
fluctuations continue to grow and so the Poisson equation 



for $> = H 



(newt) 



in the Newtonian gauge 



c K k 2 <P = 4:TrGa 2 J2Pi S 



(rest) 



(5) 



becomes dominated by the matter component, i.e. the 
dark energy is relatively smooth compared with the mat- 
ter 



Pe 6i^ « PT 4 r 



(rest) 



(6) 



where "T" denotes all other components excluding the 
dark energy. Here ck = 1 — 3K/ k 2 where K is the back- 
ground curvature. 

This condition for smoothness is not the same as set- 
ting all dark energy perturbations to zero which causes 
inconsistencies between the four scalar Einstein equa- 
tions. In particular, in the synchronous gauge, where the 
dark matter momentum also vanishes, some care must be 
taken even at ku c s 3> 1 since the dark energy momentum 
is no longer negligible in comparison Q. 

When w e = — 1, p' e / p' e w iH generally diverge leading 
to an instability in the evolution of perturbations if c 2 
is held fixed [2j. The problem arises since the change 
in the time slicing required to reach the rest or constant 
field gauge becomes infinite when the field has no kinetic 
energy. Viewed as a fluid, the problem is that the relative 
fluid velocity v e = u e /(l + w e ) becomes undefined if the 
momentum remains finite. 

If the dark energy is a composite of fields then c 2 
need not itself be fixed by fundamental properties of the 
scalars at the crossing. For example if the dark energy 



were composed of the sum of minimally coupled fields 
each with sound speed c 2 then the pressure fluctuation is 
described by 
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, (7) 



which implicitly defines c 2 as a function of the individual 
momenta. As long as no individual component crosses 
the phantom divide w ea ^ — 1, the pressure fluctuations 
are no longer singular. 

Simple two field models which cross w e = — 1 were 
constructed in 0, 3- Unfortunately, this construction 
is cumbersome for obtaining a general function w e (In a) 
constrained to match cosmological distances. 

The spirit of this construction is more broadly appli- 
cable. Models that cross the phantom divide must have 
internal degrees of freedom to ensure u e remains finite 
through the crossing. Provided they do, energy momen- 
tum conservation and the requirement that the dark en- 
ergy is smooth compared with the matter for c e kn ^> 1 
impose nearly unique constraints on their parameteriza- 
tion. We will use these requirements to construct a PPF 
description of dark energy crossing. 



III. PPF DESCRIPTION 

The PPF description of dark energy replaces the den- 
sity and momentum components with a single joint dy- 
namical variable T but retains strict conservation of en- 
ergy and momentum in its equation of motion. 

Given the conservation laws, PPF and more generally 
any minimally coupled dark energy parameterization re- 
quires two closure conditions to complete the system [U . 
The first can be taken as a condition on the anisotropic 
stress. For scalar fields this quantity vanishes for linear 
field perturbations. 

In the discussion above, the second condition was taken 
to be the relationship between pressure and density fluc- 
tuations. We saw that this choice leads to difficulties 
in parameterizing models that cross the phantom divide 
due to the appearance of singularities in the equation of 
motion for the momentum density. 

The PPF description replaces this condition on the 
pressure perturbations with a direct relationship between 
the momentum density of the dark energy and that of 
matter on large scales and a transition scale under which 
the dark energy explicitly becomes relatively smooth. 
The latter implicitly describes the momentum density 
on small scales. The strategy for choosing these relation- 
ships is to match the evolution of the metric exactly for 
scales much larger and much smaller than the transition 
scale. 

Let us start with the T variable. The conditions that 
the anisotropic stress of the dark energy vanishes and 
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FIG. 1: PPF vs scalar field calculation of the evolution of the 
potential responsible for gravitational redshifts and lensing 
($ - *)/2 for a w e = -0.7 model (flat, with Q m = 0.31 and 
h = 0.64). Curves are normalized to the initial curvature £<. 
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the Poisson equation is normal on small scales reduces 
the defining equation to (see Q Eq. 30) 



AnGa 2 . 
-t-PtIXt 

c K k z 



where 



A T = 4 re8t) = St + 3u T /k H 



(8) 



(9) 



is the density fluctuation in the zero momentum (total 
matter or comoving) gauge of the matter excluding the 
dark energy. Comparing this relationship with the Pois- 
son equation {5]) yields 



AttGcl 2 
k 2 c K 



PeSi rest) 



(10) 



The condition that the dark energy becomes smooth rel- 
ative to the matter in their respective rest gauges then 
becomes a direct requirement on the evolution of T. 

Now let us examine the second closure relation. On 
large scales, energy and momentum conservation deter- 

(T) 

mine that the curvature £ = H L in the total matter 
gauge is conserved up to order k\ in a flat universe with 
adiabatic fluctuations [l(|. The corresponding evolution 
equation for the Newtonian potentials $ and \& is closed 
by the anisotropic stress assumption [111, 03 ■ 
The Einstein equation governing £ reads 



a = i 



V T 



A-kG U e 
{aH) 2 k H lP Pe k^ 

where V T = B (T) , £ = and U e = uP 
and 



A P7 



tCkPtKT 



Pt + VT 



(11) 

in this gauge 
(12) 



FIG. 2: PPF vs scalar field calculation of the CMB anisotropy 
power spectrum for the w e = —0.7 model of Fig. [Tj and a 
two-field crossing model that approximates Wo = —1.15 and 
w a = 0.5 (flat, with Q. m = 0.26 and h = 0.74). 



with ttt as the anisotropic stress of the total matter and 

f T} 

Apr = Sp T . Since Vt = 0{knC,) we can enforce this 
condition on large scales by parameterizing a relationship 
between U e and Vt at kn <C 1 



lim U e 



H 2 



12nGp, 



■CKk 2 H V T fc , 



(13) 



where /f(lna) is a function of time only i.e. U e = 
0(kjjQ. Note that U e is the dark energy momentum 
relative to the frame defined by zero matter momentum. 
The scaling requirement is that to first order in kn, the 
dark energy and matter rest frames are the same at large 
scales. Both the single and multiple scalar field equations 
exhibit this property given that Ap e / Pe and £ are 0{k 2 H Q 
(see [lH Eq. 115 for an explicit expression). Once U e 
and its evolution are determined, 5p e follows by momen- 
tum conservation with no singularities encountered as w e 
crosses the phantom divide. 

The PPF description can be made an exact match at 
large scales to any given system of scalar fields with an 
arbitrary equation of state evolution u> e (lna) by solving 
the full equations at kn — ► and inferring for the 
evolution of all other finite k modes. However for the 
purpose of obtaining the correct evolution for the metric 
or gravitational potentials, even this is not necessary as 
long as < p e /{pT + p e )- By construction, the metric 
condition £' = 0(k 2 J ) is satisfied and the specific value 
chosen just determines the ratio of the dark energy to 
matter contributions to the metric fluctuations. Since 
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ultra large scales where the dark energy is not smooth 
are generally probed gravitationally via gravitational red- 
shifts and perhaps lensing in the future, it suffices for 
most purposes to simply take = 0. 

The final piece in the construction is to assure that the 
dark energy becomes smooth relative to the matter in- 
side a transition scale c e kn = 1 while exactly conserving 
energy and momentum locally by taking 

(l+clk 2 H )[T' + T + clk 2 H T]=S, (14) 

where 

S = [fdPT +PT) (Pe+Pe)} -|| • (15) 

This relation explicitly guarantees that T <§; Vr/kjj = 
C((f>) for crkn ^ 1. Comparison with Eq. |(5J) shows that 
this condition requires the dark energy to be smooth rel- 
ative to the matter (see Eq. . While the specifics of how 
rapidly the dark energy becomes negligible in contribut- 
ing to gravitational potentials below this scale depend 
on the specific form of Eq. (fTl"]) . the net impact on ob- 
servable quantities of this choice is small as we shall see 
below. 

The main task is to calibrate the scale of the transition, 
i.e. a relationship between cr and c e . We find that 

c r = 0.4c e (16) 

matches the evolution of scalar field models. We show 
an example with w e = —0.7 of the evolution of the quan- 
tity (<& — v P)/2 that is responsible for gravitational red- 
shifts and lensing in Fig. [TJ Metric evolution for scales 
c e kn <C 1 and c e kn 3> 1 show exact agreement between 
the PPF prescription and the direct scalar field calcula- 
tion by construction. In this model the two limits differ 
by 44% in the fractional change in the gravitational po- 
tential during the acceleration epoch. 

In Fig. [2j we compare the CMB temperature power 
spectrum in the PPF approximation to the direct scalar 
field calculation for the w e = —0.7 model and a two field 
model that approximates w e (\na) = wq + (1 — a)w a with 



wo = —1.15 and w a = 0.5 0- The latter model has w e 
evolves across the phantom divide. 

IV. DISCUSSION 

We have shown that the PPF prescription for describ- 
ing the evolution of metric perturbations in an FRW uni- 
verse is sufficiently general to encompass multiple scalar 
field models whose joint equation of state evolves across 
the phantom divide at w e = —1. This description is ac- 
curate to well below the cosmic variance limit as long as 
the transition scale to relative smoothness is comparable 
to the horizon. Moreover it is in fact exact for the metric 
evolution well above and well below the transition scale. 
As such it provides a well-controlled approximation for 
any model where the energy and momentum of the dark 
energy is separately conserved and features a transition 
of this type. 

This prescription is useful for the joint analysis of 
growth and distance measures of the dark energy, espe- 
cially those involving horizon scale perturbations like the 
integrated Sachs-Wolfe effect in the CMB. The CAMB 
Einstein-Boltzmann package has been altered to include 
PPF [lj] and a version for the dark energy has been 
made publically available [l6|]. Potential future uses in- 
clude principal component approaches to dark energy 
constraints where w e is allowed to cross the phantom 
divide multiple times (e.g. [H[). Here explicit matching 
to multiple scalar fields is cumbersome if not impossi- 
ble. The PPF prescription provides a simple but general 
approach that explicitly enforces conservation of energy 
and momentum and all of the Einstein equations remov- 
ing potential ambiguities to the meaning of a "smooth" 
dark energy component. 

WF is supported by the Initiatives in Science and En- 
gineering program at Columbia University and by DOE 
contract DE-FG02-92ER-40699; WH by DOE contract 
DE-FG02-90ER-40560, the David and Lucile Packard 
Foundation and the KICP under NSF PHY-01 14422; AL 
by an STFC Advanced Fellowship. 



arXiv:astro-ph/0407107 



W. Hu, 



[1] A. Vikman, 

[2] 
[3] 
[4] 



Phys. Rev. D71, 



Phys. 



Rev. D71, 



023515 (2005), 
047301 (2005), 

[arXiv:astro-ph/0410680] . 

R. R. C aldwell and M. Doran, P hys. Rev. D72, 043527 
(2005), [arXiv:astro-ph/0501104] . 
B. Feng, X.-L. Wang and X.-M. Zhang, 
B607, 35 (2005), 



arXiv:astro-ph/0404224 



Phys. Lett. 



[5] G. Huey, |arXiv:astrc-ph/0411102| ' 
[6] W. Hu and I. Saw icki, Phys. Rev. D76, 104043 (2007), 

|arXiv:0708.1190| . 
[7] W. Hu, P hys. Rev. D77, 103524 (2008), 

larXiv:0801.2433| . 
[8] W. Hu, Astrophys. J. 506, 485 (1998), 



[9] 



[10] 
[11] 



Mukhanov 
Lett. 85, 



and P. 

4438 



arXiv:astro-ph/9801234| . 
C. Armendariz-Picon, V 

hardt, Phys. Rev . 

[arXiv:astro-ph/0004134| . 

J. M. Bardeen, Phys. Rev. D 22, 1882 (1980). 
W. Hu and D. J. Eisenstein, Ph ys. Rev. D 59 
(1999), [arXiv:astro-ph/9809368l . 
[12] E. Bertschinger, Ast rophys. J. 648, 797 



Stein- 
(2000), 



083509 



(2006), 



arXiv:astro-ph/0604485 
[13] W. Hu, ICTP 14, 149 (2003), [arXiv:astro-ph/0402060] . 
[14] W. Fang et al, [arXiv:0808. 22081 

[15] D. Huterer and G. Starkman, Phys. Rev. Lett. 90, 



031301 (2003), [arXiv:astro-ph/0207517 
[16] http://camb.info/ppf/ 



